Home | About Me | Final Project
Geologic maps serve as a valuable tool in resource exploration. Traditional mapping is characterized by time-intensive and high-cost fieldwork, but through multispectral and hyperspectral remotely sensed images it is possible to map extensive areas instantly.
Being able to calculate areas of interest for potential hydrothermal alteration zones can be extremely useful in determining the history and available resources geologically through remote sensing.
Background
Bandpass Filters
Bandpass filters are essential for collection and analysis of satellite imagery, they allow specific ranges of electromagnetic wavelengths, or bands, to be transmitted while blocking others.
FWHM
Full Width Half Mass (FWHM) is an important concept for understanding how different bands are delineated.
The graph below visualizes the idea of FWHM. Each band in a spectral device only takes in a specific range of wavelengths, often the wavelength for each band is just denoted as the peak, but in reality each sensor is collecting a range of wavelengths. This is determined by denoting the FWHM as the area between the 50% transmittance on both sides of the curve.
Landsat 8, Hyperion, and ASTER Bands
In this plot you can see all the bands for Landsat 8, Hyperion, and ASTER. With their FWHM denoted as boxes. Its important to point out how Landsat 8 does have a band in the Short Wave Infrared (SWIR), but the FWHM is wide and this means that any SWIR data saved for this band could be attributed to any wavelength inside of the FWHM.
Location of Research
The Marysvale volcanic field in Utah is one of the largest volcanic fields in the western United States, known for its uranium and alunite deposits. The igneous rocks found here were intruded ~32-20 million years ago, due to the subduction of the Farallon plate, and extension of the Basin and Range.
Landsat
Classification of Land Cover
It is important to eliminate areas that are not exposed soil or rock to accurately map alteration zones.
Land cover classification over North America is available for download from the National Land Cover Database. 2021 classification was used as that is the most recent year available.
The Landsat imagery tested against the NLCD was collected Oct 19th, 2023.
2021 NLCD was the most recent published dataset. I reduced the dimensions to three classifications: Barren, Snow, and Vegetation
Original RGB of Landsat 8 Image
Classification Tree
I selected 50 pixels of each class that I designated as Barren, Vegeations, and Snow. These were used in a CART (Classification and Regression Trees) algorithm.
Classified Landsat Image
Confusion Matrix
Using the NLCD as a reference model, I created a confusion matrix to determine the accuracy of the created predicted raster.
Class 1: Barren
Class 2: Snow
Class 3: Vegetation
## Confusion Matrix and Statistics
##
## Reference
## Prediction 1 2 3
## 1 253987 95 280759
## 2 10001 37 607
## 3 55484 176 399263
##
## Overall Statistics
##
## Accuracy : 0.653
## 95% CI : (0.6521, 0.654)
## No Information Rate : 0.6804
## P-Value [Acc > NIR] : 1
##
## Kappa : 0.3326
##
## Mcnemar's Test P-Value : <2e-16
##
## Statistics by Class:
##
## Class: 1 Class: 2 Class: 3
## Sensitivity 0.7950 1.201e-01 0.5866
## Specificity 0.5875 9.894e-01 0.8259
## Pos Pred Value 0.4749 3.476e-03 0.8776
## Neg Pred Value 0.8593 9.997e-01 0.4842
## Prevalence 0.3193 3.079e-04 0.6804
## Detection Rate 0.2539 3.698e-05 0.3991
## Detection Prevalence 0.5346 1.064e-02 0.4547
## Balanced Accuracy 0.6913 5.548e-01 0.7063
Sensitivity (True Positive Rate): The proportion of actual positives correctly predicted.
Specificity (True Negative Rate): The proportion of actual negatives correctly predicted.
Pos Pred Value (Positive Predictive Value): The proportion of predicted positives that are actually positive.
Neg Pred Value (Negative Predictive Value): The proportion of predicted negatives that are actually negative.
Prevalence: The proportion of actual instances of each class in the dataset.
Detection Rate: The proportion of actual positive instances correctly predicted.
Detection Prevalence: The proportion of predicted positive instances.
Balanced Accuracy: An average of sensitivity and specificity, useful when classes are imbalanced.
Masked Landsat Image
Using classification as a mask to keep only ‘barren’ pixels.
Mineral Group Mapping
Through specific combinations and equations of bands it is possible to create classification rasters for varying mineral groups and assemblages. To avoid over-classification I created a threshold which I deemed as only keeping the highest potential pixels.
Iron Oxides
Iron Oxides mapping can be created using Band4/Band2.
Clay and Hydroxyl
Hydroxyl mineral mapping can be created using Band6/Band7
Ferrous
Ferrous (iron minerals) can be mapped using Band6/Band5
Sabins Ratio
A Sabins Ratio plot can be created by combining these three plots into an RGB profile.
Iron-oxide dominated areas are mapped in pink, clay and hydroxyl minerals in green and ferrous minerals are discriminated in blue. Hydrothermal alteration areas are represented by the association of green-pink or yellow zones.
ASTER
Calcite Index
Calcite minerals can be mapped using (b6/b8) * (b9/b8)
Kaolinite, Sericite, Chlorite, and Epidote Index
Kaolinite, sericite, chlorite, and epidote minerals can be mapped using (b7+b9)/b8